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Abstract 

The dark matter of our galactic halo may be constituted by elementary particles that interact 
weakly with ordinary matter (WIMPs). In spite of the very low counting rates expected for 
these dark matter particles to scatter off nuclei in a laboratory detector, such direct WIMP 
searches are possible and are experimentally carried out at present. An introduction to the 
theoretical ingredients entering the counting rates predictions, together with a short discussion 
of the major theoretical uncertainties, is here presented. 
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This is a phenomenological introduction to the detection of dark matter through its scat- 
tering in a laboratory detector. For dark matter in the form of massive quasi-stellar objects, 
like brown dwarfs, which are much bigger and much heavier than the Earth, this type of de- 
tection is quite impracticable if not undesirable. I therefore consider dark matter in the form 
of elementary particles. 

Many particles, most of which hypothetical, are at present candidates for dark matter: 
neutrinos, neutralinos, axions, etc. The methods employed in hunting for these particles are 
very different. In this short note I focus on this meeting's category of particle dark matter, viz. 
weakly interacting massive particles or WIMPs. 

WIMPs, in a broad sense, are particles with masses of the order of atomic masses or higher 
(m ^ lOGeV/c^) that interact with ordinary matter with cross sections typical of the weak 
interaction or smaller [a ^ 10~^^cm^ off a proton). The presently most popular WIMP is 
the yet-undetected neutralino, the lightest supersymmetric particle in supersymmetric models. 
Other famous WIMPs are Dirac and Majorana neutrinos, which however, thanks to the on-going 
dark matter searches complemented by accelerator results, we know not to be the dominant 
component of our galactic halo. 

A general introduction to dark matter has been given by Olive at this meeting. Direct 
detection of WIMPs was first explored by Goodman and Witten.§) General reviews are Primack 
et all-* and Smith and Lewin.i-' Engel et al^ present the nuclear physics involved. At this 
meeting Cabrera discusses experimental aspects of direct dark matter detection, while I focus 
on the theoretical aspects. 

It is worth recalling some properties of the dark halo of our galaxy. Even if recent obser- 
vations might change the details of our picture, the 1981 model by Caldwell and OstrikerB is 
good for my purposes. The Sun lies at a distance of ~ 8.5 kpc on the disk of our spiral galaxy, 
and moves around the center at a speed of ~ 220 km/s. The luminous disk extends to ~ 12 
kpc, and is surrounded by a halo of ~ 100 kpc where globular star clusters and rare subd- 
warf stars are found. Dynamical arguments suggest that the halo is filled with dark matter, 
whose local density in the vicinity of the Sun is estimated to be pdm = 0.2-0.4GeV/c^/cm^ 
= 0.7-1.4 10~^^g/cm.^ Equilibrium considerations also give the root mean square velocity of 
halo consituents to be 200-400 km/s, not much different from the escape speed from the galaxy 
(500-700 km/s). Very little is known on the mean rotation speed of the halo, and we will as- 
sume it does not rotate. All in all, there is an optimistic factor of 2 uncertainty on the density 
and velocity of halo dark matter. 

Are there WIMPs in our galactic halo? The scientific way to answer this question is to 
detect them. Signals could come from WIMP annihilation (indirect detection) or WIMP scat- 
tering (direct detection). In the former we search for rare annihilation products like neutrinos, 
antimatter or gamma-ray lines. This is reviewed by Bergstrom at this meeting. In the latter, 
the basic philosophy is to build a target, wait and count. 

The WIMP scattering rate per target nucleus is the product of the WIMP flux 0^ and of 
the WIMP-nucleus cross section a^^i. For an order of magnitude estimate we take the effective 
WIMP-nucleon coupling constant to be Fermi's constant Gp = 2.3016 10~^^^^c^ cm/GeV, 
which sets the scale of weak interactions. We distinguish two cases: (i) the WIMP couples 
to nucleon spin, a^^i ~ G^nl/ff ^ 10~^'^cm^ ; and (ii) the WIMP couples to nucleon number, 
a^i ^ Gpfi'^A^/h'^ ^ 10~^°cm^. Here /ij = m^nii/ {m^ + mi) is the reduced WIMP-nucleus mass, 
and Ai is the atomic number of the target (^ 80 in the numerical examples) . The WIMP flux is 
0^ = vp^jm^ ^ 10'^cm~^s~^/(m^c^/GeV), for a WIMP density ^ 10~^^ g/cm^ and a typical 
WIMP velocity v ~ 300 km/s. The resulting scattering rates, taking ^ 100 GeV/c^, are of 
the order of ^ 1/kg-day for spin-coupled WIMPs and of ^ lO^'/kg-day for WIMPs coupled to 
nucleon number. These rates are quite small compared with normal radioactivity background. 



Therefore the common denominator of direct experimental searches of WIMPs is a fight against 
background. 

For this we get help from characteristic signatures that we do not expect for the background. 
For example, while the Earth revolves around the Sun, the mean speed of the WIMP "wind" 
varies periodically with an amplitude of 60 km/s. This leads to a ~ 10% seasonal modulation 
in the detection rate, with a maximum in June and a minimum in December.^* As another 
example, the direction of the WIMP "wind" does not coincides with the Earth rotation axis, 
so the detection rate might present a diurnal modulation due to the diffusion of WIMPs while 
they cross the EarthB* (this however occurs for quite high cross sections). A final example of 
background discrimination is that the WIMP signal is directional, simply because most WIMPs 
come from the direction of the solar motion.S* 



WIMP-nucleus scattering 

Since the relative speed v ^ 300km/s ~ 10~^c, the process can be treated non-relativistically. 
The center of mass momentum is given in terms of the reduced WIMP-nucleus mass as = /XjV 
and is ^ AiMeV/c since fii < rrii. The corresponding de Broglie wavelength is ^ 200fm/Aj, 
and can be smaller than the size of heavy target nuclei, in which case nuclear form factors are 
important. In the laboratory frame, the nucleus recoils with momentum q = 2/;; sin(^cm/2) and 
energy u = g^/2mj. Here ^cm is the center-of-mass scattering angle. The 4-momentum transfer 
is very small, ^ 10~^ GeV^/c^ (compare with a typical deep inelastic ^ 1 GeV^/c^). 
The differential scattering rate per unit recoil energy and unit target mass is formally 



The sum is over the nuclear isotopes in the target, Tj(g^) is the scattering matrix element at 
momentum transfer squared q"^ = 2miZ/, and fi is the mass fraction of isotope i. A sum over 
final and average over initial polarizations is understood in |rj(g^)p. The factor 

V^{q) = r (2) 

Jq/2fj.i V 

with units of inverse velocity, incorporates the x velocity distribution f^iv). For a Maxwellian 
distribution with velocity dispersion frms, seen by an observer moving at speed vq, 



2vo 



(3) 



with Vq = q/2fii. For standard halo parameters, r]i{q) is approximately exponential in the 
deposited energy u. The previously-mentioned modulations enter the rate through rii{q). 
The scattering matrix element T{q^) can be written as the Fourier transform 

T(g2) = j (f|\/(f)|i)e'^^/Mf (4) 
of a non-relativistic WIMP-nucleus potential 

nr)= E {G: + Gy^a^)5{r-r^). (5) 



pointlikc 
nuclcons 
n— p,n 



The constants and G'2 are effective four-fermion coupling constants for nucleon-WIMP inter- 
actions, and are analogous to Fermi's constant Gp. represents scalaiQ or spin-independent 
interactions, axial^ or spin-dependent interactions. Both terms are coherent in the quantum- 
mechanical sense when g-Rnucieus ^ ^, when the nucleus can be treated as pointlike and 
T(g^) can be taken as T(0). At larger q, which can occur with heavy target nuclei, both terms 
are incoherent. Nuclear form factors F{q^), conventionally defined by T{q^) = T{0)F{q'^), 
should then be introduced. The scalar and spin form factors are in general different, reflecting 
the difference in the mass and spin distributions inside the nucleus. 

The task of a theoretician is to provide a theoretical estimate of T(g^) starting from a 
particle-physics model. We accomplish this by stages, successively finding the WIMP-quark, 
the WIMP-nucleon and the WIMP-nucleus effective lagrangians. Step 1, finding the effective 
WIMP-quark lagrangian at small g^, is analogous to going from the Standard Model to four- 
fermion interactions. Step 2 requires knowledge of the quark content of the nucleon, i.e. the 
contributions of different quarks to the nucleon mass and spin. Step 3 needs a nuclear model 
to describe how protons and neutrons are distributed in a nucleus. 

This procedure is now illustrated for a Dirac neutrino and for a Majorana particle, an 
example of which is the neutralino. 



Dirac neutrino 

Step 1: a Dirac neutrino u interacts with a quark q through the diagram in Fig. la. At <^ m|, 
the Z propagator reduces to ig'^'^/m'^, and the four-fermion amplitude reads 

^/2GF ^{vu - a^75)7^z/ q{y^ - aq75)7''q, (6) 

with Vy = ay = \, a^ = T^^ and fq = — 2eq sin^ 6w Here sin^ 6w — 0.23 and Cq and T^q are 
the electric charge and the third component of the weak isospin of quark q. For a non-relativistic 
neutrino, only the time component of the vector current and the space components of the axial 
current survive. The first is spin-independent {u-'^qu oc u^p) and the second spin-dependent 
(1/775!/ oc u^au). 

Step 2 for the vector part 

V2Gf VyVq z/7^i/ q7''q, (7) 

because of vector current conservation, simply amounts to summing T^q and of the constituent 
quarks. For protons and neutrons one obtains respectively 

GP = ^{l-Asin'OwW (8) 
G: = -^vy. (9) 

The interaction is mainly with the neutrons since 1 — 4 sin^ 9w ~ 0. 
Step 2 for the axial part 

V2Gf a^aq 2/7^751/ q7^75q, (10) 
leads to the four-fermion coupling constants 

Gl = V2GFau{anAu + adAd + asAs) , (11) 
Gl = V2GFau{a^Ad + adAu + a^As) . (12) 



^Associated to scalar and axial vectors under 3d rotations. 




Figure 1: Examples of WIMP-quark scattering. 



Aqs^. It can be 



1.2573 ±0.0028 



Here Aq is the fraction of the proton spin carried by quark q, |(p|q7^75q|p) ■ 
obtainedtl'' from data on neutron and hyperon /3-decay, which give Au — Ad 
and Au + Ad — 2As = 0.59 ± 0.03, respectively. The contribution of the strange quark is 
As = in the naive quark model, As = —0.11 ± 0.03 ± ■ • • from deep inelastic data, and 
As = —0.15 ± 0.09 from elastic z/p z/p data. 

Step 3 for the spin-independent part introduces the nuclear mass form factor -Finass(5'^), and 
results in 

massVi ) 



(13) 



where N {Z) is the number of neutrons (protons) in the nucleus. Neutron scattering off nu- 
clei suggests that F^assil'^) — F^.ra.i.'f) ^ the electromagnetic form factor. The electric charge 
distribution is well-described by a Fermi or Woods-Saxon form,0^ whose Fourier transform is 
indistinguishable from the convenient analytic expression0) 



qR 



(14) 



The electromagnetic radius R and the surface thickness s can be obtained by fitting electron 
scattering dataJl^P or can be roughly approximated by i? ~ A^f^ fm and s ~ 1 fm.o' -Fe.m.(5'^) 
presents diffraction zeros when the modified Bessel function ji{qR) = 0, the first of which 
occurs at qR ~ 4.2. In electron scattering, these diffraction zeros are filled in, because due 
to the long-range Coulomb attraction the electron wave function is distorted from a simple 
plane wave and the form factor is not simply the Fourier transform of the charge density. The 
short-range nature of WIMP-nucleus interactions make us expect no wave function distortion, 
and diffraction zeros remainP] The first diffraction zero is important in assessing bounds from 
some present-day detectors.ll3) 

Step 3 for the spin- dependent part requires the expectation values of the total spin of protons 
{Sp) and neutrons {S^) separately. At g = 0, 
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(15) 



where J is the nuclear spin. Even-even nuclei, with even numbers of protons and of neutrons, 
do not have spin, and for them Ta{0) = 0. For even-odd nuclei with J 7^ 0, a nuclear model is 

^The first corrections are at a level of 10~^ and come from neglected higher powers of the incoming WIMP 
velocity. 



The single-particle shell mo 



needed to estimate (Sp) and {SnL For instance, '^Ge is an odd-neutron nucleus with J 
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0.50, 



(5p) = 0; 



(16) 



the odd-group modeljil^ in which the odd-nucleon spin is related to the nuclear magnetic 
moment fi and gyromagnetic factors g^p, gives 



9n-9n 



0.23, 



a more sophisticated interacting shell mode© gives 

(Sn) = 0.468, (^p) = 0.011. 



:i7) 



;i8) 



The proton might have a small but non-zero contribution to the cross section, which might 
change the relative merits of different nuclei for dark matter searches. 

At g 7^ 0, nuclear spin form factors are needed. The neutron and proton contributions 
differ, and at present only complex calculationslll' 0* for specific nuclei provide an estimate of 
the isoscalar and isovector spin form factors F^^^^^q"^) and F^pi^^q"^), in terms of which 
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(GI + G:) (Sp + S^) F^.{q') + (GI-G:) (Sp-S^) F'{q 



(19) 



The results of these calculations can be conveniently resumed by the approximate expressions 
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with parameters given in the following table for selected nuclei: 



J 


^'max/keV 


(5p) 




ro/fm 


ri/fm 


c/fm 


valid for 
u/kcV< 


28Si 1 

39k ^ 

2 


540 
216 
100 
145 


0.011 
-0.0019 

0.3430 
-0.184 


0.468 
0.133 
0.269 
0.054 


1.971 
1.302 
1.378 
1.746 


2.146 
1.548 
1.600 
1.847 


-0.246 
-0.320 
0.196 
0.371 


55 
145 

'^max 



Major ana fermion 

A Majorana fermion is a spin-| particle that coincides with its antiparticle. It carries no 
conserved quantum number. It has neither vector nor tensor currents. Of the remaining 
pseudoscalar, scalar and axial currents, only the last two have a non-vanishing non-relativistic 
hmit, spin-independent the first {xx ^ X^x) and spin-dependent the second (x775X ^ X^^x)- 

Axial currents may arise from exchange of a Z boson as in fig. lb, and the analysis is then 
analogous to that in the previous section, with the obvious replacement of a,y with a^. 

Scalar currents originate from exchange of a scalar particle ip, e.g. as in Fig. Ic. At small 
g^, the ip propagator reduces to —i/rn? and the four-fermion amplitude reads 



(21) 
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Figure 2: Scattering rate versus mass for neutralinos: (a) phenomenological approach,t3' (b) grand- 
unified approach^^ 



For a nucleon n = p, n one then obtains 



m 



(22) 



For example, in the case of the neutrahno with exchange of the lightest supersymmetric Higgs 
boson, the sum over quarks is explicitly 
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cos a, _ _ sma, _ , 

(muUU + mcCC + m^tt) -{m^dd + mgSS + mbbb) 



sin P 



cos (3 



(23) 



The scalar quark content of the nucleon (n|qq|?7,) can be extracted from data with the help of 
chiral perturbation theory, vr-nucleon scattering and heavy quark expansion.0) The result is 



(muuu) ~ (mddd) ~ 30MeV/c^, 
(mcCc) = (mbbb) = (mttt) = ^ | rrip 



(m^ss) ~ 60-120 MeV/c^ 
E ("^qqq) I - 60MeV/c^ 

q=u,d,s 



(24) 
(25) 



The strange quark contribution is uncertain by a factor of 2. Step 3 is analogous to the Dirac 
neutrino case, and leads to eq. (|13|) with four-fermion couplings given by (p2|). 



Neutralino 

Supersymmetry and the neutralino have been presented by Jungman at this conference. The 
neutralino has both spin-dependent and spin-independent interactions with nuclei, the former 
mediated by Z boson and squarks, the latter by Higgs bosons and squarks. The general for- 
malism of the preceding sections can be used. In the limit of heavy squarks q^, the effective 
four-fermion constants are given by 

GP~G-= ^ (qq)|- ^ ^/.XX^/»qq ^l^ ^Lq,xq^i?q.Xq j ^ ^^6) 
q=u,d,s,c,b,t \ h=Hi,H2 k=l "^qfe / 



GP= 51 Aq (" ^^xxf qq + 1 ^ ^^q^-xq + ^kxq \ ^ g° = gp(Au ^ Ad). (27) 

q=u,d,s V "^z 8 ^^-^ m-^ y 

Expressions for the elementary vertices gijk can be found in ref. 18. 

Predictions in supersymmetric models suffer from the presence of many unknown parame- 
ters. Two extreme attitudes are a phenomenological approach in which what is not excluded 
is allowed, and a grand-unified approach in which coupling constants and masses are unified 
at some high energy scale. Fig. 2 shows examples of calculated event rates in ^^Ge, each point 
representing a choice of model parameters: "predictions" may well span 10 orders of magnitude 
in a phenomenological approached and 2 orders of magnitude in a more restricted scenario.li2P 



Underabundant dark matter relics 

Given a particle-physics model, the relic density of a species, a WIMP x particular, is 
a calculable and definite quantity. Often it happens that the computed relic density is 
(much) smaller than the dark matter density in the Universe. For this reason, some authors 
simply neglect this case. But even if these WlMPs constitute only a fraction of the dark matter, 
they generally have quite high scattering cross sections off nuclei, because of an approximate 
inverse proportionality of the x I'elic density and the x^nucleus cross section. However, the 
scattering rate also includes the x halo density p^. It is reasonable that is only a fraction 
of the local dark matter density pdm, but which precise fraction it is depends on the model for 
galaxy formation. If both the main and the x components of dark matter are cold, we expect 
them to behave similarly under gravitation, so that the halo fraction might be equal to the 
universal fraction Q^/Qdm- Unfortunately, ^dm is poorly known: it can range from ^ 0.01 for 
dark matter associated with galactic halos to ~ 1 for a smooth universal component. In fig. 3, 
the suppression of scattering rates due to rescaling of the neutralino halo density by a universal 
fraction with Qdm^^ = 0.025 is apparent to the left of the dashed line. This suppression must 
be included for consistency when setting bounds on particle-physics models. 
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Figure 3: Scattering rate versus relic density for neutralinos (from ref. 18). 
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